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We discuss properties of the quadrupole collective excitation of the deformed neutron- 
rich nucleus 38 Mg within the framework of quasi-particle random phase approximation 
(QRPA). We first solve the coupled-channels equations to obtain the single-particle levels, 
and construct the ground state by treating the pairing correlations in the BCS approx- 
imation. We then solve the QRPA equation using the response function formalism, by 
including the continuum spectra with the box dicscretization method. We show that the 
collectivity of the gamma vibration (the lowest K n = 2 + mode) is significantly enhanced 
if protons and neutrons have different deformations. We also discuss an attempt towards 



full continuum QRPA calculations for deformed nuclei. 



1. Introduction 



Without an exception, stable nuclei with a neutron number iV=20 or 28, as well as their 
neigbouring nuclei, have a spherical intrinsic shape. Heavy neutron-rich nuclei, however, 
often exhibit deformed characters even around the shell closure. Typical examples in- 
clude 32 Mg [ [T] , 34 Mg [ 12 Ej , and 44 S [ E] . This phenomenon is often interpreted as a 
disappearance of spherical magic shells [ E] • 

In this contribution, we investigate the excitation structures of such deformed neutron- 
rich nuclei within the framework of quasi-particle random phase approximation (QRPA). 
For nuclei far from the f3 stability line, it has been well recognised that both the pairing 
interaction and the couplings to the continuum spectra play an important role. The most 
robust and consistent framework to deal with these effects on the ground state properties 
is the Hartree-Fock-Bogoliubov (HFB) method [ [7J [H] - The continuum QRPA framework 
based on the HFB approximation has also been developed recently [ EH] • A problem 
with the HFB method, however, is that it can be computationally very demanding. For 
this reason, the application of the continuum QRPA on top of the HFB ground state has 
been limited only to spherical systems so far [Ellin]- In this paper, we avoid this difficulty 
by using the BCS approximation. This is an approximation to the HFB method, where 
one takes only time-reversed pairs in the Bogoliubov transformation. With this approxi- 
mation, the particle-hole and the particle-particle channels are decoupled completely (for 
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a density-independent pairing force at least) both in the ground state and in the QRPA 
calculations. Notice that the BCS method has been shown a reasonable approximation 
to the HFB, within the same treatment for the continuum states, even for nuclei close to 
the neutron drip line [ |Hl HI] • 

Another difficulty related to a deformed QRPA calculation is that the configuration 
space can be quite large. In the past, QRPA calculations for a deformed nucleus have 
been done either with a simple separable interaction [ or with a drastic truncation 
of the configuration space [ ■ in this paper, we present a formalism which uses the 
response function in the coordinate space [E], which is most suitable to a zero-range 
interaction, eliminating these limitations. Together with the BCS approximation for the 
pairing interaction, we will demonstrate below that a deformed QRPA calculation now 
becomes feasible within a reasonable computation time. 



2. Formalism 

2.1. Coupled-channels approach to single-particle spectra 

Let us consider a single-particle motion in a deformed mean-field potential, 

V(r) = V ccnt (r) + V(V ls (r)) ■ (-iV x tr), (1) 
with 

V cent (r) = Vo{r) - RP^p- Y2o(r ) + (v c (r) + 3{Z ~ 1} *PY M {*)) (1 - r z ), (2) 

dV (r) 

VUr) = VM-RsoP^^Y^r). (3) 

Here, (3 is the deformation parameter. For simplicity, we take into account only the linear 
order of (3 in the potential. We have also assumed an axial symmetry, where both the 
parity ir and the spin projection K onto the z axis are conserved. 

In order to obtain the single-particle levels and corresponding wave functions, we solve 
the coupled-channels equations for both negative and positive energy states [US]. This 
method has an advantage that the correct scattering boundary conditions can be easily 
implemented for the positive energy solutions, that is important in a full continuum 
calculation for neutron-rich nuclei. In this paper, for the sake of simplicity, we instead 
discretize the positive energy spectra within the same framework by putting the nucleus 
in a large spherical box, as was recently demonstrated in Ref. [ [Tj)] , and leave the full 
continuum treatment SIS db future work (see sect. 4). 

The coupled-channels method has sometimes a difficulty in finding all solutions in a 
deformed potential, especially when there are two states close to each other in energy. One 
can easily miss a few eigenvalues in a search routine. We have overcome this difficulty 
by introducing the generalized wave function node for a multi-channel system, which was 
originally proposed by Johnson We have checked the efficiency of this method by 

applying it to a spherical system. We have also computed the Nilsson diagram of 238 U 
and have confirmed that all the levels obtained were smoothly connected to the spherical 
levels. Evidently, the coupled-channels method with the generalized wave function node 
efficiently works even for heavy nuclei, where there are lots of intruder levels. 
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Using the single-particle levels obtained in this way, we then solve the BCS gap equation. 
As we emphasized in the introduction, we use the BCS approximation rather than the 
HFB approximation. The applicability of the BCS approximation for neutron-rich nuclei 
has been well tested in Ref. [ |HJ EJ. In this paper, we particularly use the density- 
dependent zero-range force for the pairing interaction. We will specify the parameters in 
the next section. 

2.2. Response function formalism for deformed QRPA 

The elements of the QRPA matrices have been given in the textbook of Rowe [ E] : 

Aap^s = (E a + Ep)5 al3a s + (a8\v ph \pi) (u a vp + v a u p ) (u y v s + v y u s ) 

+ (aP\v pp \j5) (uaUpu^us + v a VpVjV S ), (4) 

B af3r/ s = (aj\v ph \^5) (u a Vf3 + v a up){u 7 v s + v^u s ) 

-(a(3\vpp\j5) (u a upv 7 v s + v a vpu 7 u s ), (5) 

where E a is the quasi-particle energy and a is the time reversed state of a. In writing 
down these equations, we have used the identity, (a5\v\/3~f) = — {a^\v\^8), which follows 
from the property of the time-reversal operator. Notice that the particle-particle matrix 
elements vanish in the BCS approximation unless (3 = a and 5 = 7 in eqs. (0J) and (0). 
On the other hand, the particle-hole contribution is finite only if (3 7^ a and 5 7^ 7 due to 
the BCS phase convention, Va = —v a and Ua = u a . Therefore, the particle-particle and 
the particle-hole channels are decoupled completely in the BCS approximation, not only 
in the ground state but also in the QRPA calculations. 

Retaining only the particle-hole channel, one can write down the QRPA equations in 
the following form: 

= (E a + E p -u)X aP + J drdr'D a0 (r)v ph (r,r')T(r'), (6) 
= (E a + E p + u)Y aP + J drdr'D* a - p {r)v ph {r,r')T{r'). (7) 

Here, we have explicitly expressed the matrix elements of the particle-hole interaction in 
the integral form. Using the single-particle wave function ^ Q (r), D a p(r) and T(r) are 
defined by 

D a p(r) = (u a v/3 + v a ui3)(f>* a (r)(f>p(r), (8) 
T(r) = 52(D;- s (r)X lS + D y - s (r)Y jS ), (9) 

respectively. From these equations, one obtains the generalised (Q)RPA dispersion rela- 
tion 



= Jdr' 5(r-r')-f dr" U (r,r") v ph {r" ,r') T{r') = J dr' W{r,r')T{r'), (10) 

where the unperturbed response function is given by 

_ / D*Jr)D ad (r') D aB (r)D* -Jr') \ 

II (r, r") = -£( = fj - + = • (U) 

^ \E a + Ep - 10 - %r\ E a + E p + u + it] j 
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The response of the system to an external field can be calculated using the QRPA 
response function defined by [ IT4] 



n(r, r') = J dr" W 



(12) 



For axially symmetric wave functions (fi a , the response function can be decomposed into 
multipoles as 



n(r,0 = J]^ — U K (rL,r'L') Y LK (r)Yl, K (r). 



(13) 



L.L' K 

We assume a zero range particle-hole interaction of the form 



v p h{r, r 



.X:even 



(14) 



Then, the QRPA response function (JT2jl is obtained by solving the Bethe-Salpeter equation 
for each K, 



Il K (rL,r'L') =IlP(rL,r'L') + J2 / dr" [n (0) i 

L" 



K 



(rL,r"L")U K (r"L",r'L') 



where 



U^v] K (rL,r'L') = Y,T, U P( rL ^' L >^ r ')( Y L''K\Yx \Y LIK ). 



(15) 



(16) 



A L" 



The strength function S for an external field V ex t(v) = f^Yix^) is then obtained in 
the usual way as, 



S{u>) 



--Im f drdr'V; xt (r)Tl(r,r')V cxt (r') 

7T J 

If 1 

Im / r 2 drr' 2 dr' — Ilx(rL,r'L)f(r)f(r') 

7r J rr' 



(17) 
(18) 



where Im denotes the imaginary part. 

A typical dimensionality of the QRPA response function goes as follows. If one 
includes the radial mesh up to r=15 fm with 0.5 fm step, there are 30 points for the r 
integrals. For the interaction given by eq. (jl4j) . only those L with the same parity are 
coupled to each other, and generally the convergence is rather quick. Therefore, for the 
quadrupole response for example, one would include only L = 0, 2 and 4. By treating 
proton and neutron separately, the total dimension of the response function is 30 (for r) 
x3 (for L) x2 (for t z ) = 180. This is much smaller than the number of 2 quasi-particle 
configurations, which can be more than several thousands for a deformed system. 

3. Quadrupole collectivity of 38 Mg 

We now apply the formalism to the quadrupole response of the 38 Mg nucleus. We 
assume a Woods-Saxon form for the mean- field potential Vo(r) and V so (r) in eqs. (j2J) and 
(J2J) with parameters given in Ref. [I19j. We take the deformation parameters (3 from the 



5 



results of HFB calculations with a Skyrme interaction by Terasaki et al. [I2DJ- For 38 Mg, 
the deformation parameters are 0.33 and 0.28 for protons and neutrons, respectively. 
Notice that the deformation parameters are different for protons and neutrons by an 
amount of 0.05. We use the same pairing interaction as in Ref. [ I2UJ for our BCS 
calculations. This is a density-dependent zero-range interaction with an energy cut-off 
at 5 MeV above the Fermi energy. Although a smooth energy cut-off was employed in 
Ref. [1201, we use a sharp energy cut-off for simplicity. With these parameters, we obtain 
A p = —24.4 MeV and A n = —0.56 MeV for the proton and the neutron Fermi energies, 
A p = 0.0 MeV and A n = 1.50 MeV for the proton and the neutron average pairing 
gaps. For the residual interaction v p h, we use the SIII parameter set [ l2*T] of the Skyrme 
interaction with the Landau-Migdal approximation for the momentum dependent terms [ 
ITS] . We include the single-particle levels up to 50 MeV in the calculations. The continuum 
states are discretized with the box discretization method with the box size of 15 fm. We 
have checked that the results are not altered by more than 10% even if we use a larger 
box size. For single-particle states above the pairing cut-off energy, we simply take v 2 = 
for the BCS occupation probabilities. 

Since our calculations are not self-consistent in the sense that the residual interaction 
is not related to the mean field potential, we renormalize the residual interaction as / ■ v p h 
so that the Goldstone modes appear at zero excitation energy. For a spherical system, 
there is only one Goldstone mode for the center-of-mass motion, and one can use a single 
value for the renormalization factor in order to define the effective residual interaction and 
study several multipole responses. In contrast, for a deformed system, the center-of-mass 
motion splits into two excitation modes according to the K quantum number (K 71 = 0™ 
and 1~). Moreover, there is another Goldstone mode associated with the rotation motion 
of the whole nucleus. This motion carries K* = 1 + . We find that the renormalization 
factors required to put these modes at zero energy are significantly different from each 
other for the configuration space mentioned above. The renormalization factors are found 
to be 1.00075 for K" = 0", 0.95585 for K n = 1", and 1.0300 for K* = 1+ Since we are 
interested in only the quadrupole motion in this paper, our prescription is to use fi+ for 
the other K quantum numbers, i.e., K w = + and 2 + . Note that the Goldstone mode 
for the number conservation does not appear in our calculation, since we include only the 
particle-hole channel in QRPA. 

Figure 1 shows the isoscaler (IS) quadrupole strength distribution of the 38 Mg nucleus 
(the lower panel) in comparison with that of the 24 Mg nucleus (the upper panel). We set 
r] = 0.5 MeV in eq. (jllj) in order to smear out the results for the presentation purpose. 
The dotted, dashed, and solid lines are for K w —0 + , 1 + , and 2 + modes, respectively. Since 
we use a density-dependent pairing interaction, there is a small coupling between the 
particle-hole and the particle-particle channels for the K n = + mode (notice that a time- 
reversed pair can form only a K n = + state) even in the BCS approximation, when the 
residual interaction is evaluated as the second derivative of the energy functional [IU] ITU]. 
We simply neglect this effect in this work. The K 7T =1 + and 2 + modes are free from this 
approximation. 

In the figure, one observes the splitting of the giant quadrupole resonances (GQR) [ 22J 
in the energy region 10 - 20 MeV both in 24 Mg and in 38 Mg. In contrast to the stable 
nucleus 24 Mg, we notice that the IS GQR is more fragmented in the neutron-rich nucleus 
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38 Mg due to a large neutron contribution. We have also computed the neutron response, 
and found that the K 7 * = 2 + mode in the energy region 3-10 MeV is almost a pure 
neutron mode. Another interesting result is that the collectivity of the lowest K n = 2 + 
mode is increased considerably in 38 Mg. This is due both to the pairing effect [E] and 
to the different deformation between proton and neutron, as we discuss below. 

In Ref. [ 2(J|, Terasaki et al. have suggested that a low-lying quadrupole isovector mode 
may appear when the protons and the neutrons have different deformation parameters. 
In neutron-rich nuclei, the isovector and the isoscalar modes are coupled to each other [ 
and thus this effect would also enhance the low- lying isoscalar quadrupole strength. 
In order to confirm this, we repeat the calculations by setting the proton and the neutron 
deformation parameters to be identical. Figure 2 shows the low-lying part of the isoscalar 
quadrupole response of 38 Mg for the K 11 = + (the top panel), K n = 1 + (the middle panel), 
and K n = 2 + (the bottom panel). The dotted and the dashed curves are for the case 
where the proton and the neutron deformations are set to be the smaller value 0.28 and 
the larger value 0.33, respectively. The solid line denotes the original calculation shown in 
fig. 1, where the proton and the neutron have different deformations. For K 77 = + and 
1 + , those three lines are very similar to each other, while, in marked contrast, the lowest 
K n = 2 + excitation (the gamma vibration) clearly shows a different behaviour. When the 
deformation is different between proton and neutron, the excitation energy of this stae 
becomes significantly smaller and at the same time the strength is increased considerably. 
The effect of a larger proton deformation parameter simply changes the dotted curve to 
the dashed curve. There is an additional effect originating from the different proton and 
neutron deformations, that enhances the collectivity of the gamma vibration. 

4. Concluding remarks 

At present, a continuum QRPA calculation based on the HFB method is still too costly 
to apply to deformed nuclei. The only possible way to overcome this difficulty is to em- 
ploy the BCS approximation for the pairing correlations, where the particle-hole and the 
particle-particle channels are decoupled. Contrary to a general belief, the BCS approxi- 
mation for the description of ground state properties has been shown adequate even for 
neutron-rich nuclei [ |Hl HI] • We presume that this is the case also for the description of 
excited states, and have developed the QRPA method based on the BCS approximation. 
Particularly, we have presented the response function method for deformed QRPA, which 
substantially reduces the computational effort to solve QRPA compared with the conven- 
tional method which diagonalises the non-hermitian QRPA matrix. We have applied the 
method to quadrupole responses of the deformed neutron-rich nucleus 38 Mg. We have 
shown that the collectivity of the gamma vibration is increased when the protons and the 
neutrons deform in a different way from each other as the mean-field calculation predicts. 

In this paper, we have discretized the continuum states by putting the nucleus in a 
large spherical box. For a more consistent description, one would need to treat the 
continuum spectra exactly. This applies both to the BCS and to the QRPA calculations. 
In the continuum BCS calculation for the ground state, the single-particle resonance 
states play an essential role [E]. Recently, we have shown that the coupled-channels 
approach with the eigenchannel basis offers an efficient way to identify a deformed single- 
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particle resonance state through the so called eigenphase sum We have also derived 
a factorization formula for the energy and the radial dependences of a scattering wave 
function near an isolated resonance state [US]. This formula will be useful when one 
computes the matrix elements of pairing interaction over a wide energy region including 
the continuum spectra. As for a continuum QRPA calculation, it is well known that 
the exact treatment of the continuum couplings is possible if one uses the single-particle 
Green function in the coordinate space representation [121] • An extension of this approach 
to a deformed system has been done in the field of atomic and molecular physics [ |2H] > 
which again uses the coupled-channels method with the eigenchannel basis (see also Ref. [ 
)• We have also recently extended the Green function approach to spherical open-shell 
nuclei within the BCS approximation [ 14. We plan to combine these two extensions of 
the Green function approach in order to calculate the continuum response of superfluid 
deformed nuclei built on top of the continuum BCS ground state. We will report the 
results elsewhere. 
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Figure 1. The QRPA responses of the 
24 Mg nucleus (the upper panel) and of 
the 38 Mg nucleus (the lower panel) to an 
isoscalar quadrupole field. The dotted, 
dashed, and solid lines are for the K n = 
+ , 1 + , and 2 + modes, respectively. 



Figure 2. The effect of different proton 
and neutron deformations on the isoscalar 
quadrupole response of the 38 Mg nucleus. 
The dotted and the dashed lines are ob- 
tained by setting the proton and the neu- 
tron deformations to be the same as in- 
dicated in the inset, while the solid line 
is the result when the proton deformation 
is larger than the neutron deformation by 
0.05. The top, middle, and bottom panels 
are for the K w = + , f + , and 2 + modes, 
respectively. 



